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Abstract 

In this report we summarize the results of a three year research program on high 
pressure vapor transport (HPVT) of compound semiconductors -funded -by-NASA Grant- 
N-AGW-2865-Most of our work focused onto pnictides, in particular ZnGeP 2 , as a model 
system. Access to single crystals of well controlled composition of this material is desired 
for advancing the understanding and control of its point defect chemistry in the context of 
remote, real-time sensing of trace impurities, e.g. greenhouse gases, in the atmosphere by 
ZnGeP 2 optical, parametric oscillators (OPOs). The report is presented in 9 parts 
addressing 

\ 1. Background information concerning the motivations lor our research 

\ ' \ 

\ 2. Experimental arrangements for high pressure vapor transport 

\ \ 

\ 3. Identification of transport reactions and thermochemical calculations 

\ \ 

\ 4. Phase equilibria and transport conditions \ 

5. Estimation of physico-chemical parameters for the modeling of HPVT 

6. Modeling of flow under the conditions of HPVT 

7. Characterization of crystal properties - \ 

8. Point defect chemistry and physical properties 

9. Interactions and presentations 

A comprehensive assessment of the current status of. the art and science and ZnGeP2 


crystal growth and characterization was carried out as part of an international workshop 

held in Vesser, Germany, August 21-26, 1995 ^ A brief summary of the results of this 

workshop is included as an addendum to this report (attachment A-l). 
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1. Introduction 

The transport of transparent crystals of ZnGeP 2 in a temperature gradient using a liquid 
source was discovered by Xing et al. in 1989 2 and suggested superior optical properties, 
that is, lower residual absorption in the transparency range below the fundamental 
absorption edge than for melt-grown single crystals of the same material. This is of 
interest in remote detection, for example, in the context of continuous emission 
monitoring (CEM) devices, relying on continuously tunable infrared lasers, such as, 
ZnGeP 2 optical parametric oscillators (OPOs). Minimum absorption is essential in this 
application for optimizing the laser damage threshold and the efficiency of frequency 
conversion. The bandgap (~2.2 eV) and infrared absorption edge (~0.1 eV) of ZnGeP 2 
makes it possible to cover the wavelength regime between 2-13 ptm where many 
fundamental rotational-vibrational molecular states are found. At present, the fabrication 
of ZnGeP 2 OPOs is based on bulk single crystals grown from the melt for which, in spite 
of 25 years of research and development, the absorption behavior is still poorly 
understood. Our research concerning the HPVT growth of ZnGeP 2 crystals and their 
characterization aims at improving this situation by the generation of standards of well 
defined point defect chemistry - based on the far more sensitive control of the solidus 
composition by the control of partial pressures of the constituents during vapor growth 
than possible under the conditions of growth from the melt. 

2. Experimental conditions of HPVT 

Figure 1(a) - (c) show schematically the experimental arrangements used thus far for the 

study of HPVT growth in our laboratory. The arrangement used by Xing et al. 2 is shown 

in fig. 1(a) and involved transport in a shallow vertical temperature gradient from molten 
ZnGeP 2 that is contained in a pyrolytic BN boat to the wall of the fused silica enclosure. 
Subsequent studies employed a zone melting geometry that allowed growth from a 
moving molten source zone established in a bar of ZnGeP 2 thus minimizing changes in 
the composition of the transported crystals due to depletion of the melt in the course of 
the growth process. A typical crystal obtained in this way is shown in figure 2 3 . It is 
transparent and of bright red collor. On its surface it exhibits a set of macrosteps due to 
step bunching in the course of the growth process. A third arrangement is similar to that 
described by Scholz 4 and employs a ring shaped source heated by glow bars to 
temperature T s with transport to a window on the bottom of the fused silica enclosure 
held at temperature T c 5 (see attachment A- 2). The regions above and below the glow bar 
position are kept at constant temperature T, by means of two isothermal furnace liners 
filled with sodium vapor and the window is cooled by a jet of nitrogen gas to establish 


well controlled thermal conditions. A booster heater is inserted in the top well of the 
enclosure to permit fine tuning of the temperature T p . 

3. Identification of the transport reaction and thermochemical calculations 

The partial pressures of P 2 /P 4 vapor and zinc over ZnGeP 2 have been established through 
preceding work( 6 M 8 ) and are shown as a function of temperature in table 1. To our 
knowledge the nature of the Ge vapor species participating in the transport reaction has 
not been addressed in the literature. The vapor pressure over pure germanium in the 
temperature range of ZnGeP 2 HPVT is well known, that is, < 2x1 O' 8 bar at T < 1300K. 
Since the partial pressure of atomic Ge vapor over molten ZnGeP 2 is expected to be less 
than that over pure Ge at the same temperature the observed transport rates, typically 
5xl0 -6 mole/m 2 s, cannot be explained in terms of atomic germanium vapor playing a 
significant role in the transport reaction. Therefore, molecular Ge containing species must 
be present in the vapor phase. A possible source for the generation of such volatile 
molecular species is residual oxygen. Thermochemical calculations using the Chemsage 
code 9 were carried out, for the oxygen partial pressure expected on the basis of (in 
decreasing order of importance): (i) oxygen diffusion through the fused silica wall of the 
growth capsule, (ii) outgasing of residual water from the silica walls and (iii) reaction of 
the vapor atmosphere with the inner wall of the capsule. They predict a partial pressure of 
GeO(g) of 10' 12 <pGeO<10~ 9 in the temperature range 1100K<T<1300K 10 (see 
attachment A-3 for details). However, cutting off the supply of oxygen by diffusion 
through the fused silica wall by the use of nitrogen jacketed fused silica ampules, and 
reducing out-gasing as a source of oxygen by a careful vaccum bake of the growth 
capsules, does not significantly reduce the transport rate. 

Thus the observed efficient transport must be explained by an alternative transport 
mechanism, that is, molecular Ge x P y must be present in the vapor phase over molten 
ZnGeP 2 - This is corroborated by optical absorption spectroscopy 11 . Figure 3 shows a 
schematic representation of the experimental arrangement used by us for this study. 
Fused silica cuvettes ranging from 50 to 100 mm length and 10 mm diameter, containing 
red phosphorus as well as ZnGeP 2 , ZnP 2 and GeP, respectively, were heated in a furnace 
- equipped with isothermal furnace liners to establish a uniform temperature over the 
length of the cuvettes. The absorption spectra for the vapor atmosphere above these 
materials were measured in the range 180 < X < 1pm. The details of the measurements 
and resulting spectra for the four materials investigated at X < 350nm are given in the 
attached reprint A-4 of ref. 11. A comparison of the spectra for the three compounds 
under study shows spectral lines associated with vibronic transitions in the vapor phase 



above ZnGeP2 that cannot be explained by the absorption lines for either P2, P4, Zn, their 
oxides or GeO. 

Table 2 shows the position of the lines associated with Ge x P y in comparison to 
spectral lines associated with P2. Although the resolution of our spectrometer is not 
sufficient to allow an unequivocal identification of the stoichiometry of the Ge x P y vapor 
species, based on the limited number of known Ge-P compounds, we identify it 
tentatively as GeP, which is presently being further investigated by mass spectrometric 
analysis. Thus the reactions 

Zn(g) + GeP(g) + l/2P2(g) <=> ZnGeP2(s) 
and 

Zn(g) + GeP(g) + l/4P4(g) <=> ZnGeP2(s) 
are chosen to describe the overall transport. 

4. Phase equilibria and transport conditions 

Figure 4 shows the Gibbs phase triangle for the system Zn-Ge-P. Of the Ge-P 
compounds shown on the Ge-P binary only GeP, and possibly GeP 2, exist at 1 bar 
pressure. GeP3 and GePs can only be prepared at very high pressure ( 65 kbar) 12 . On the 
Zn-P binary two compounds, Zn3P2 and ZnP2 exist at lbar pressure. Using the 
thermochemical data of table 3 , the predominance field for ZnGeP2 has been calculated 
and is shown in fig. 5 for two temperatures that are relevant for the execution of HPVT 
crystal growth. Note that the limit at very high phosphorus pressure is not well defined by 
the calculation. However within the range of phosphorus pressures used in our ex- 
periments the width of the predominance field ZnGeP2 clearly decreases with decreasing 
temperature and the limits of the predominance field shift toward smaller Zn pressures. 
Also, the predominance fields shown in fig. 5 allow ZnP2, Zn3P2 and ZnGeP2 or Zn3P2, 
Ge and ZnGeP2 to coexist in equilibrium. 

Table 4 shows the results of transport experiments employing ZnGeP2 charges plus 
various additions, e.g., O2, phosphorus, ZnP2, and Zn3P2, respectively. Also, the 
temperatures of the source and of the growing crystal was varied within the ranges 130 IK 
< T s < 1315 K and 1248 K< T c < 1280 K. Shallow temperature gradients (~ 4 K/cm) and 
short transport distances are required for single phase transport by HPVT. Thus the 
vertical transport observed in ref. 1 corresponds to close to ideal conditions. Also, a set of 
experiments was conducted metering well controlled amounts of HC 1 gas into the growth 
capsule to enhance the Ge-concentration in the vapor phase. Halide transport of ZnGeP2 



has been utilized by others 13 for homoepitaxial growth on ZnGeP2 substrates of (100), 
(001) and (112) orientation according to the transport equation 


2Zn(g) + P2(g) + GeCl2(g) <=> ZnGeP2(s) + ZnCl2- 

Single phase ZnGeP2 HC1 transport is attained for a typical composition of 6.27xl0' 7 
mole/cm 3 of zinc, 3.254xl0' 5 mole/ cm 3 phosphorus and 1.3xl0‘ 5 mole/cm 3 HC1. 


5. Estimation of physico-chemical parameters for the modeling of HPVT 

The relative importance of convective and diffusive transport in HPVT is captured by the 
solutal Peclet number that is defined as 



where v c , l c , and D represent the linear flow velocity, a characteristic length scale, in this 
case the diameter of the vortices, and diffusivity. The flow is characterized by the 
Reynolds number 
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where p is the vapor density and r\ is the kinematic viscosity. For large values of the 
Peclet number boundary layers form and the transport becomes diffusion limited. Under 
conditions of laminar flow, that is, Re < 2100, the width of the boundary layer is 
proportional to 1/VRe. 

Based on a preceding evaluation of the second virial coefficient for phosphorus 
vapor 14 , in the pressure range considered for HPVT of ZnGeP2 deviations from ideal 
behavior are small so that the kinetic theory of gases provides valid order of magnitude 
estimates of D and r|. At 1300K, the vapor atmosphere over ZnGeP2 consists of 
approximately 2.77 bar of P2 and .38 bar of P4 corresponding to a mass density p =5.17 
kgnr 3 . The sum of the radii of P2 and P4 rp2-P4 = 5.2xl0' 10 m resulting in a = 8.5xl0~ 19 
m 2 and A, = kT/V2op = 4.75xl0' 8 m. With reduced molecular mass m r = 6.9x10 ~ 26 kg the 
mean molecular speed c av = (8kT/7tmr) 1/2 =814 ms -1 , that is, r) ~ (nm r Ax av NL)/3V « 
1.56xl0" 5 kgnrV 1 and DP2/P4 - A,c av /3 ~ 1.29xl0 -5 m 2 s _1 . The heat capacity at constant 
pressure C p = 606 Jkg^K* 1 Kt = (nA,c av C v )/3 = 0.227 JK^s^nr 1 . Also, under the 



conditions of our experiments, l c ~ 10' 2 m and v c ~ 5xl0' 2 ms _1 , so that Pe(ideal gas) = 39 
and Re(ideal gas) = 172. 

Using corrections for molecular interactions in the computations of the viscosity, 
thermal conductivity, and diffusivity 15, 16 r| = 2.69xl0' 5 kgnrV 1 (269 micropoise), K t = 
2.30xl0’ 2 JK'W and D P2 p 4 (1300K) = 2.26xl0’ 5 mV 1 . Thus Pe d = 22.1 and Re = 96, 
that is, for the case of ZnGeP 2 , both convective and diffusive transport have to be 
considered, but the flow can be expected to be laminar, which simplifies the 
computations. In order to take into account the temperature variations in the HPVT 
enclosure we consider next the Rayleigh number 

Ra = p 2 gcttCpAT 
rpct 

with regimes of flow 

laminar: Ra<10 8 

transitional: 10 8 < Ra <10 10 

turbulent: Ra > 10 10 

Under the conditions of our ZnGeP 2 experiments, that is AT = 80 K and the volume 
expansion coefficient oc t ~ nR/p = 2.64xl0* 5 m 3 K‘ 1 . Thus Ra = 1.12xl0 5 so that the 
previous prediction regarding the flow based on the Reynolds number holds. Although 
for ZnGeP 2 the flow thus is predicted to be laminar more detailed numerical calculations 
are needed for its quantitative characterization under the conditions of HPVT 
experiments with the specific thermal boundary conditions for the geometries represented 
in Fig. 1, taking into account the thermal properties of the fused silica enclosures (K t = 
1.381 Js- 1 m^K-l). 


6 . Modeling of the flow and real-time process monitoring 

The mathematical model for transport processes in chemical vapor deposition are 
described by conservation of mass 
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Based on the conservation of mass it can be shown that 



dV = 0 


J[|p + div(pu)] 

for any domain Q t of the gas phase. Thus the integrand must vanish in accord with the 
conservation of mass equation given above. For incompressible fluids this reduces to 



div(u) = 0 . All buoyancy effects are included in the model through the gravitational term 
in the Navier-Stokes (momentum) equations that are based on Stokes' hypothesis relating 
the shear viscosity to the bulk viscosity so that the viscosity can be represented by a 
single parameter r) 16 . Also, body force enters into the transport model through the 
momentum equations and the separation of gas species due to different size and mass is 
modeled by the terms 3/3^[pDAB«TfAfB()logT73 £] in the transport equation, where 0.7 is 
the thermal diffusion constant. In the Boussinesq approximation 17 , the density is assumed 
to be a constant in all equations except for the existence of a buoyancy force in the 
presence of a gravitational field due the density variations related to thermal expansion, 
representing the body force term pg as p[l - Px(T-T re f)g, wherep is a constant, Pt is the 
thermal volume expansion coefficient and T re f is the reference temperature. 

The coupled set of non-linear differential equations thus defined was solved 
numerically for a range of parameters which include: reactor geometry, operating 
pressure, gravitational vector, and the parameters evaluated in section 5 , that is, Kt, q, p, 
andC v . Two fluid dynamics packages, FLUENT, a finite difference code, and FIDAP, a 
finite element code, were licensed for the simulation study. Due to the ability to model 
more complex geometries the FIDAP package was found to be more versatile for our 
work and most of the computations were based on this code. Initially the flow dynamics 
of a homogeneous gas inside a closed ampule of fused silica was modeled numerically 
for axi-symmetrical conditions using quadratic-triangular meshes. The numerical 
simulations were performed to study the flow dynamics and temperature distribution 
inside the fused silica ampule and to illustrate the feasibility of an optimal reactor design 
study. Several parameter studies were caried out . The problem governed by the coupled 
set of partial nonlinear differential equations was simulated for a range of parameters 
with a view to maximizing transport and uniformity of deposition. These parameters 
include: reactor geometry, aspect ratio, operating pressure, gravitational vector, thickness 
of the fused silica window and boundary temperature profiles. The parameter studies 
indicate thet the best configuration has the substrate on the bottom and that flow rate 
increases with increase in pressure ' . In addition, conditions for optimum convective 
mixing at the source and transport to the seed crystal in a closed ampule high pressure 
vapor transport were obtained, and led to changes in the design, including the location of 
the heating element 5 . 

In addition to the parameter studies, a more rigorous approach to the optimal 
design and control has been taken. In ref. 20 the experimental design problem was 
formulated as a shape optimization problem with respect to the geometry of the reactor 
and a boundary temperature control problem. The material and shape derivatives of 



solutions to the Boussinesq approximation were derived. Optimality conditions and a 
numerical optimization method based on the augmented Lagrangian method were 
developed for boundary control of the Boussinesq flow. The regularity of solutions has 
also been studied. In ref. 21; well-posedness and necessary optimality conditions for 
optimal control problems for homogeneous incompressible gas flow inside a HPVT 
ampule by boundary temperature control was studies. A second order method based on an 
augmented Lagrangian formulaton was applied to solve the optimal control problem for 
the steady flow. Our numerical calculations indicate the effectiveness of temperature 
control through a portion of the boundary for improving the vertical transport flow inside 
the ampule. 

The modeling was then further advanced to permit the investigation of non-axi- 
symmetric geometries where 3-dimensional modeling must be applied 22 (see attached 
preprint A-4b for details). The latter is mandatory in the context of real-time optical 
monitoring of the growth process for which we introduced a new method that is highly 
sensitive with regard to both the initial phase of nucleation and late stages of three 
dimensional growth 23 . This method, p-polarized reflectance spectroscopy (PRS), is based 
on the changes in the reflectance of the overgrowth of a substrate with regard to a beam 
of p-polarized light that impinges onto the surface at the Brewster angle of the substrate 
as illustrated in Fig. 7 for the overgrowth of a seed crystal plate (e.g., a wafer of silicon) - 
located on the cooled fused silica window of the growth ampule depicted in Fig. 1(c) - by 
a crystal layer (e.g., ZnGeP 2 ). The reflected beam is detected by a photodiode, processed 
through a phase sensitive amplifier and read into a computer. 

The complex reflectivity r of a three-layer stack composed of materials labeled a, c and s 
for ambient, crystal and seed, respectively, with interfaces labeled ac and cs, is given by 
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where e a , £c, and e s refer to the dielectric functions at the energy of the incident laser 
beam of the ambient, the crystal and the seed assuming that reflections at the seed crystal 
interface with the cooled fused silica window can be neglected. Due to constructive and 
destructive interference between the partial waves reflected at the crystal/seed and am- 
bient/crystal interfaces, upon initiation of growth the reflected intensity R p = rr* oscillates 
between maxima and minima the distance of which in time is uniquely related to the 
phase thickness of the growing crystal layer that enters into the periodic exponential 
function in the equation for r (see Figure 8 for a simulation of the growth of various 
materials on Si). Thus the growth rate should be accessible in real-time if the problems 
related to beam stirring in a dense vapor atmosphere can be controlled. In addition to 
providing information regarding the growth rate, a fine structure superimposed on the 
oscillations due to interference is observed for pulsed vapor deposition processes. The 
analysis of this fine structure results in detailed information on the kinetics of crystal 
growth that is, on the properties of the surface reaction layer that connects the growing 
film to the nutrient vapor phase. Initial simulations of this four-layer stack have been 
carried out for the growth of GaP on Si 24 (see Attachment A- 10). 

In order to assess this problem for the conditiond of HPVT of ZnGeP 2 numerical 
calculations of the density variations and the related gradients in the dielectric function of 
the ambient phase, the vapor, were performed using non-uniform quadrilateral elements 
for the representation of the 3-D flow in the growth capsule that no longer supports axi- 
symmetric flow. This is due to the alteration of the capsule shown in Fig. 1(c) by the 
provision of two additional optical windows for the implementation of PRS, that is, an 
entrance port for the laser beam to strike the substrate located on the bottom of the 
capsule at its Brewster angle and an exit port for the reflected beam that is intercepted by 
a detector. Under the conditions of close-space vapor transport of ZnGeP 2 from top to 
bottom at the earth's gravity and 10" 2 g, which represents an upper limit for the g-jitter 
expected for the space shuttle during astronaut activities 22 , our computations show that 
changes in the ampule orientation with regard to the gravity vector affect the temporal 
gradients by less 16% at lg and <1% at 10 _2 g. Thus only moderate reductions of the 
gravity vector are needed for this type of experimentation and astronaut activities do not 
impose a significant limitation to their successful excecution. 

An assessment of the problems associated with beam stirring in the context of 
real-time optical process monitoring was made on the basis of combining Snell's law 



and the Dale-Gladstones law 


nismcpi = n2Sintp2 


ni =hpi + 1 


with Dale-Gladstone constant h = 2.26xl(H m 3 kg' 1 to compute the deviation of the angle 
from normal incidence at the entrance window along a path inside the growth capsule - 
divided into 26 segments ending at the point of reflection on the substrate wafer. The 
result is shown in Fig. 9 and indicates that beam bending in the relatively steep 
temperature gradient close to the substrate surface may affect the sensitivity and accuracy 
of the method. However, calculations - based on a reliable parameter set - can provide for 
corrections that restore high sensitivity and accuracy by appropriate compensation of the 
entrance angle. Also, the problem of noise due to random time-dependent fluctuations of 
the vapor density can be largely eliminated by phase-sensitive detection so that with 
sufficient computational support PRS would be a reliable method for remote real-time 
sensing of crystal growth processes on earth and in space for a wide range of vapor 
pressures. 

7. Characterization of crystal properties 

Figure 10 shows steady state photoluminescence (PL) spectra for (a) ZnGeP 2 
platelets grown by the HPVT method and (b) ZnGeP 2 plates cut from a bulk single 
crystal grown by the gradient freezing method with and without annealing and etching. 
All exhibit in the as-grown state broad infrared emission peaked at 1.3 eV 25 . The PL 
spectra for a crystal grown by HPVT from a nearly stoichiometric melt in a shallow 
gradient, in addition to deep luminescence, exhibit near edge emission. For bulk single 
crystals grown by the gradient freezing method, near edge emission is absent. However, a 
high energy shoulder is observed in the dominant deep emission extending to ~1.7 eV. 
Subband excitation extending to 1.6 eV are also observed for these crystals in the 
photoconductivity (PC) spectra. Annealing in vacuum at 500°C for 400 hrs results in a 
substantial increase of the 1.7 eV emission relative to the 1.3 eV emission. Also, a 
significant reduction in the near infrared absorption is observed upon annealing. 

Figure 1 1 shows the time dependence of the PL, which cannot be described by 
exponential functions, that is, the PL intensity has a hyperbolic dependence on time 26 . 
Thus the PL spectra are dominated by donor-acceptor transitions that are associated with 
time constants depending on the relative spacing of the donors and acceptors in the 
lattice. 



Figure 12 shows an interpretation of all optical data known todate the position of 
energy eigenvalues and subbands in the bandgap of ZnGeP 2 that are associated with 
residual donor and acceptor states. For example, the broad PL band, peaked at 1.3 eV, is 
attributed to transitions between filled acceptor (Al) and unfilled donor (Dl) states that 
are associated with subbands having a center distance of 1.3 eV. Explaining the sub- 
bandgap structure in the PC spectra by transitions between filled Al states and the 
conduction band edge (CBE), the highest filled Al state is associated with an energy 
eigenvalue ~1.6 eV below the CBE. The region of residual absorption within the 
transparency range of ZnGeP 2 extends from the absorption edge to well below 1.6 eV, 
due to transitions between filled acceptor and empty donor states that do not contribute to 
the PC spectra because of the high effective mass in the donor subbands. 

There are at least three effects that must be considered in the interpretation of the 
annealing studies: (i) The reduction of quenched-in disorder, reducing the concentration 
of point defects, (ii) The shift of the FL position associated with this annihilation of 
defects, altering the population of the remaining defect states, (iii) Shifts in the FL 
position caused by the generation of defects at the crystal surface upon annealing. The 
observation of additional structure in the PL and the depth dependence of the 
luminescence features after annealing clearly establishes the latter and requires the 
consideration of additional donor-acceptor energy subband transitions. 

8. Point defect chemistry and physical properties 

For the above mentioned thermal conditions and partial pressures, at the onset of 

the HPVT experiment, the supersaturation at the substrate location corresponds to zinc 
and phosphorus supersaturated growth conditions. However, upon prolonged HPVT 
growth without melt replenishment, the gradual depletion of the source of zinc causes 
conversion to Zn-deficient growth conditions that may result in the formation of a Ge- 
rich second phase on the surface of the ZnGeP 2 crystals in extreme cases. The changes in 
the point defect chemistry accompanying the transitions from maximum to minimum 
zinc activity in the solid are reflected in the optical properties of the ZnGeP 2 crystals. 
Thus the gradual transitions from the border of the homogeneity range about ZnGeP 2 
corresponding to maximum zinc activity to the border corresponding to maximum Ge 
activity during growth confirms our initial premise that HPVT permits the control of 
stoichiometry over the entire existence range of ZnGeP 2 - 

Point defects that may be generated upon deviation of the ternary compound com- 
position from stoichiometric ZnGeP 2 are Zn , Ge and P vacancies - denoted by V zn , V a 
and V p - and various antisite defects, e.g., Ge on Zn sites denoted by Ge m . Note that 
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rich second phase on the surface of the ZnGeP 2 crystals in extreme cases. The changes in 
the point defect chemistry accompanying the transitions from maximum to minimum 
zinc activity in the solid are reflected in the optical properties of the ZnGeP 2 crystals. 
Thus the gradual transitions from the border of the homogeneity range about ZnGeP 2 
corresponding to maximum zinc activity to the border corresponding to maximum Ge 
activity during growth confirms our initial premise that HPVT permits the control of 
stoichiometry over the entire existence range of ZnGeP 2 - 

Point defects that may be generated upon deviation of the ternary compound com- 
position from stoichiometric ZnGeP 2 are Zn , Ge and P vacancies - denoted by V Zn , V Ge 
and V P - and various antisite defects, e.g., Ge on Zn sites denoted by Ge^. Note that 



specific deviations from stoichiometry within the homogeneity range about the exact 
composition ZnGeP 2 does not establish the dominance of a specific native point defect, 
since combinations of other defects may create the same deviations and interactions with 
extrinsic point defects must be considered. 

However, certain conclusions are possible by the interpretation of the optical 
characterization results in conjunction with the results of electrical measurements and 
ESR/ENDOR investigations. For example, ZnGeP 2 and ZnSi x Ge!. x P 2 epilayers grown 
by organometallic chemical vapor deposition exhibit a wide range of solid solutions 
corresponfing to the replacement of a large fraction of the zinc atoms by group IV atoms 
in metastable equilibrium depending on the zinc to germanium precursor flow rate ratio 27 ' 
29 . This is only possible by massive antisite defect formation on both the zinc and 
phosphorus sublattices. Also, regardless of the as-grown composition, annealing in 
vacuum depletes the near surface region with regard to the most volatile constituents, that 
is, phosphorus and zinc. By the mass action law, the formation of vacancies on the zinc 
sublattice reduces the concentration of Znce acceptors and enhances the concentration of 
Gezn donors according to the reactions [1]: V Zn + Znce Zn x and [2]: V Zn + Ge x 
Gezn- The introduction of V P donors at the surface is associated with the formation of 
and Gep acceptors according to the reaction [3]: V P + Ge* -» Gep, and possibly, but not 
necessarily, of Znp acceptors. Thus the luminescence at 1.6 and 1.7 eV observed in 
vaccum annealed bulk crystals (see attachment A- 12) may be explained tentatively by 
transitions from V P donors to V Zn and Gep acceptors, respectively. The overall FL shift 
associated with the generation of V P and V Zn depends on the values of the equilibrium 
constants for reactions [1] - [3], which are not known at present. The high p-type 
conductivity observed upon annealing in zinc vapor, has been explained by Rud as an 
effect of Zn P antisite defect formation 32 . 

A more definitive control of the stoichiometry of ZnGeP 2 is obtained by growth 
under the conditions of HPVT. The thermochemical calculations of the temperature 
dependence of the existence range of ZnGeP 2 show that HPVT of ZnGeP 2 from a molten 
source ties at the lower substrate temperature into solidus compositions close to the 
boundary between the ZnGeP 2 and ZnP 2 of the ZnGeP 2 existence range 11 . The n-type 
conductivity of ZnGeP 2 crystals grown under such conditions suggests the formation of 
relatively shallow donor states that move the FL toward the conduction band edge. 
Simultaneously emission at 1.8 eV is observed that is neither seen in melt-grown crystals 
nor for growth from zinc-depleted vapor sources. The high phosphorus pressure used in 
the HPVT experiments rules out a higher V P concentration than observed under the 
conditions of melt growth. Thus the donor responsible for the 1.8 eV emission for HPVT 



growth from a zinc supersaturated vapor phase must be associated with a different point 
defect that does not form under for HPVT growth from zinc-depleted melt compositions. 
Possible choices are the P m and P Ge antisite defects of which only the P Ge donor explains 
the decrease of the intensity of the 1.8 eV emission with decreasing zinc activity in the 
solid. Therefore, we associate the relatively shallow donor that is formed upon HPVT 
growth from a phosphorus and zinc supersaturated vapor phase with the P Ge antisite 
defect. 

Further information is gained from the analysis of optical brightening upon annealing of 
bulk ZnGeP 2 crystals grown from the melt in conjunction with the PL results. The broad 
deep luminescence in ZnGeP 2 crystals observed for crystals grown from the melt and 
zinc-depleted vapor sources is due to close-spaced D-A pairs , corresponding to the broad 
donor and acceptor subbands. Accepting the interpretation of the ENDOR measurements 
of ref. 30, the deep acceptor is associated with Va,. By equation 2 it is in turn associated 
with the Ge Zn antisite defect, so that the Gezn-Vz,, close-spaced D-A pairs represent a 
possible choice. The alternative close-pairing Vp-Vzn has been proposed in ref. 31 on the 
basis of the photo-induced EPR study. Unfortunately the photo-induced EPR emloying 
HeNe laser radiation probes only a thin surface region of the ZnGeP 2 crystal, which is 
expected to be depleted of Zn and P. Thus the finding of both V P and by photo- 
induced EPR, though being in agreement with other evidence refererring to the near 
surface region, does not necessarily pertain to the optical brightening upon annealing, 
which is a bulk effect. We prefer the Gezh donor since it provides for the simultaneous 
elimination of V Zn and Gezn sites upon annealing according to the reaction [4]: Gezn — > 
Ge x + V zn and the mass action law, which may be further enhanced by a simultaneous 
reduction of cation disorder on the Ge sublattice, that is, Zn Ge according to [5] Zn Ge + V zn 
-» Zn x . The associated FL motion toward the conduction band edge that results in the 
filling of residual deep Gezn donor states by electrons further enhances the optical 
brightening. Of course, the simultaneous elimination of P- and Zn-vacancies upon 
annealing would be also possible, but there exists no obvious linkage in this case. 
Additional insights are likely to emerge from the results of irradiation studies that permit 
the controlled formation of defects and their characterization. Work on this aspect is 
being continued in the PFs laboratory in collaboration with industrial partners. 


9. Interactions and presentations 

The execution of this program greatly benefitted from interactions with Dr. Warren 
Ruderman and Dr. Ilya Tseybak of INRAD, INC at Northvale, NJ and Dr. Sebastian 



Fiechter at the Hahn-Meitner Institut Berlin, Germany. Two graduate students Mr. Gary 
Wood and Mr. Hal Castleberry made valuable contributions to our research. Mr. Wood 
completed a master thesis in Materials Science and Engineering in 1994 and Mr. 
Castleberry is expected to complete a PhD thesis in Materials Science and Engineering 
on HPVT of ZnGeP 2 in the spring of 1996. The results of our work on were reported by 
Dr. Fiechter in a paper entitled “High pressure vapor transport of binary and ternary 
compound semiconductors” at the 6 th International Symposium on Experimental Methods 
for Microgravity Materials Science, held at San Francisco in 1993. Prof. Bachmann 
presented an invited paper entitled “Heteroepitaxy of wide-bandgap ternary semi- 
conductors at the 9 th International Conference on Ternary and Multinary compounds at 
Yokohama, Japan in 1993. Prof Tran presented a paper entitled “Modeling of flow 
dynamics and its impact on the optimal reactor design problem” at the SIAM meeting on 
the Identification and Control in Systems Governed by Partial Differential in 1993 and 
an invited paper entitled Gas dynamics and transport in a high-pressure reactor under 
microgravity conditions” at the 2 nd IEEE Mediterranean Symposium on New Directions 
in Control and Automation, June, 1994 in Chania, Greece. Prof. Ito presented an invited 
paper entitled “Optimal control of thermally coupled Navier-Stokes equations” at the 
IMA meeting in 1994. Two additional papers on different aspects of “High pressure 
vapor transport of ZnGeP 2 crystals” were presented by Dr. Fiechter and Prof. Scroggs at 
the 7 th International Symposium on Experimental Methods for Microgravity Materials 
Science, held at Las Vegas in 1994. On August 25, 1995 a full day of a workshop brought 
together experts on ZnGeP 2 crystal growth and characterization from the US and Europe 
to critically assess the present status of our knowledge in this field (see attached work- 
shop summary (attachment A-l). 
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Table 1 

Partial Pressures of Vapor Constituents Over ZnGeP 2 
based on a Chemsage calculation 




Table 2 

Positions of Absorption Heads 

(absorption lined due to GeP x + absorption lines of P 2 for reference) 


Absorption, lines 
of P 2 in nni 


203.6 

205 

205.5 

206.9 

207.49 

208.83 

209.43 

210.81 

212.26 

212.86 

214.3 
215.0 

216.4 

218.64 

220.14 


226.161 


232.65 

236.76 

240.99 


Absorption lines 
of Ge in nm 
HI: 19073, 194.1, 195.4, 
197.0, 198.0, 200.0. 
201.5 I 


205.35 I 

207.58 P 

209.8 I,P 

211.07 P 

212.02 I 

212.98 P 
213.62 
214.26 P 

215.94 II 

216.8 P 
217.78 I, II 
218.74 P 
220.34 I,P 

222.59 II 

223.55 I 
226.11 I 

227.08 I 
228.37 II 

231.59 I 

232.55 
234.16 1,11 
236.73 P 

238.98 I 
240.91 p,n 

241.56 

243.2 
244.51 1 
246.77 H 



Table 3 

Thermochemical Data 

(dimensions: standard heats of formation AH° f , heats of fusion AH^ U and heats of 
evaporation AH ev in kcal/mole; standard entropies AS° f and molar heat 
capacities c p in cal/mole*Kelvin) 


Element/ AHf° 

Compound AHf U AH e 

S f ° 

C p = A+B- 
A 

10“ 3 -T+C 

B 

.• 10 5 -T~ 2 

c 

P (red) 0.00 

5.45 

4 . 051 

3.559 


P 4 (white) 16.68 

39.28 

18.28 

15.12 

- 

P (1) 17.31 

- 

25.168 

- 

T>3 1 7 K 

P 4 (g) 30.77 

66.89 

19.562 

0.162 

-3.213 

P 2 (g) 40.68 

52.11 

8 . 675 

0.191 

-0.994 

Zn (s) 0.00 

9.943 

5.096 

2.782 

0.13 

Zn (1) 1.75 

- 

7.495 

- 

- 

Zn (g) 31.15 

38.451 

4.965 

- 

- 

ZnCl 2 (s) -99.13 

26.62 

14.290 

8 . 994 

- 

ZnCl 2 (l) 2.474 

- 

2.588 

- 

- 

ZnCl 2 (g) -63.839 

66.185 

14.74 

-1.03 

- 

Ge (s) 0.00 

Ge (1) 8.82 

7 .425 

5.577 

6.596 

0 . 931 

-0.25 

T>12 1 0 K 

Ge (g) 89.34 

40.1 

7.038 

-0.873 

0.76 

GeH 4 (g) 21.68 

51.89 

14 . 939 

5.276 

-5.166 

GeP(s) -0.13 

14 . 6 

10.84 

2.7 

-1.25 

GeP(g) 28.4 

44 .43 

12.34 

- 

- 

GeO(s) -50.3 

11 . 9 

9.54 3 

3.497 

- 

GeO(g) -10.99 

53.49 

8 . 92 

0.0239 

-1.3685 

Ge 2 0 2 (g) -112 . 

75 . 1 

20. 

- 

- 

Ge 30 3 (g) -212 . 

99.3 

32. 

- 

- 

GeCl(g) 174.36 

58.72 

10.314 

-0.532 

-1 . 18 

GeCl 2 (g) -40.84 

70 . 63 

13.818 

0.052 

-0.88 

ZnP 2 -24.15 

14.4 

17.01 

3.997 

0.65 

Zn 3 P 2 -38.65 

36.2 

30.15 

6.225 

-0.36 

ZnGeP 2 -26.3 

28.22 

24.65 

3.5 

-4.0 



Chemical Transport Experiments in the Parallel Processing Furnace 



Z 113 P 2 , (2. 9mg) ZnjP 3 :Ge (5%) 

P (50mg) Ge/GeP (2%) 










DM“T 


r(V | [ 


1 


Monochromator 



-y 



Motor -Controller 


Keithley 602 
Preamplifier 


Lock-in Amplifier 
EG & G 


Power 

Supply ORIEL 


ISA HR 640 


Iris " 






Quartz Cuvettes 
// 





Xe-Lamp 

I 

ORIEL 7340 


Figure 3. Schematic representation of the arrangement used in optical absorption 


measurements 
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Figure 5. Predominance field for ZnGeP 2 at two temperatures 11 . 
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Figure 6. Calculated isotherms (solid lines A 1229K, B 1242K, C 1254K, D 1267K, 

E 1279K, F 1292K, G, 1304K, H 1317K, 1 1329K, J 1342K) and flow (small 
arrows) for HPVT transport of ZnGeP 2 in a reactor of modified Scholz geometry, 
(a) 0.1 g, (b) 0.001 g. . 
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„ Figure 7. Schematic representation of the implementation of PRS under the conditions 
of HPVT in a reactor of modified Scholz geometry. 
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• angle of incidence is the Brewster angle of the Si-substrate 

(dielectric function %= 6328 A) = (15.25 , 0.17) -> <Pb = 75.637°) 



film thickness (A) 

Figure 8. Simulation of the interference oscillations in the PRS intensity for the 

heteroepitaxial overgrowth of silicon for a selected dielectric constants of 
the epilayers 
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Figure 9. (a) Isotherms and flow velocity vectors for a close-space vapor transport 

s 

ampule incorporating ports for p-polarized reflectance spectroscopy; 
(b) Deviations of the beam angle inside the growth ampule from normal 
incidence at the entrance window due to beam stirring in the vapor 
atmosphere 
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Figure 10. Photoluminescence spectra of ZnGeP 2 : (a) Crystals grown by HPVT, (b) 
Crystals grown from the melt, top trace as-grown, middle trace annealed, 
bottom trace etched 
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igure 12. Interpretation of the observed optical behavior of ZnGeP 2 associated with 

various point defect related sub-bands and energy eigenvalues in the bandgap 






